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lattice  parameters  (translations) 

lattice  parameters  (angles) 

reciprocal  lattice  parameters  (translations) 

reciprocal  lattice  parameters  (angles) 

volume  of  unit  cell 

arbitrary  directions  in  the  lattice 

the  point  on  the  projection  representing  direction  0 

Miller  indices  of  direction  0 

Miller  indices  of  direction  0' 

lattice  repeat  distances  for  directions  0 and  0' 

values  of  the  parameters  ti.  r and  it  which  give  local  maxima  in  the  value 
of  D 

an  integer 

lattice  repeat  distances  defining  the  various  ranges  within  which  D may  fall 

Miller  indices  of  a direction  which  has  a repeat  distance  Dn 

orthogonal  directions  in  the  lattice 

Miller  indices  of  direction  01 

Miller  indices  of  direction  02 

Miller  indices  of  direction  03 

centre  of  reference  sphere  and  of  projection 

radius  of  reference  sphere 

angles  made  by  direction  0 with  directions  01.02.03  respectively 
polar  coordinates  of  point  P on  the  projection 


The  Miller  indices  un.ni.wn  arc  distinct  from  nO.rO.wD.  i/l.H.n  I etc.  No  confusion  should  arise 
since  the  former  arc  used  only  in  the  context  of  subroutine  LIMITS. 


1.  INTRODUCTION 


The  stereographic  projection  offers  a convenient  method  for  mapping  certain  types  of 
crystallographic  information  in  two  dimensions  and  for  determining  angular  relationships 
between  various  features  in  a crystal.  Particularly  useful  are  “standard  projections"  of  a crystal. 
These  are  stereographic  projections  showing  the  relative  positions  of  certain  lattice  plane  normals 
and/or  crystal  directions  for  given  orientations  of  the  crystal,  especially  orientations  of  high 
symmetry.  Use  of  the  stereographic  projection  in  crystal  physics  is  described  by  Barrett  [1], 

2.  PROGRAMS  FOR  STEREOGRAPHIC  PROJECTIONS 

Many  projections  have  been  published  over  the  years  but,  except  for  the  simplest  crystal 
classes,  a projection  produced  for  one  situation  is  unlikely  to  be  of  use  in  another,  and  the  need 
frequently  arises  for  projection  along  some  unusual  direction  in  the  crystal,  for  a projection  of  a 
crystal  which  has  non-cubic  symmetry  and  for  which  a projection  is  not  already  available,  or 
for  a projection  in  some  particular  format. 

Two  computer  programs  have  been  written  to  fulfil  these  needs.  One  program  determines 
the  coordinates  of  the  points  on  a standard  projection  in  a format  specified  by  input  parameters. 
The  other,  using  output  from  the  first,  plots  the  projection  on  a graph-plotting  terminal. 

Similar  programs  have  been  written  and  published  in  the  past  [2, 3,4.5]  and  the  present  pro- 
grams perform  basically  the  same  function  as  these.  However,  the  options  written  into  the 
present  programs  provide  a combination  of  features  which  the  author  has  found  convenient  in 
deformation,  texture  ai  I electron  diffraction/microscopy  studies.  The  programs  have  been 
developed  and  tested  over  a period  of  many  years  and  include  several  novel  features. 

The  programs  are  written  in  FORTRAN-IV  as  implemented  on  a DECsystem-10.  The 
coding  has  been  largely  optimised  in  the  most  deeply  nested  sections  of  the  programs.  The 
coordinates  program.  SPC.  should  be  readily  adaptable  to  any  system  with  a FORTRAN-IV 
compiler  since  an  attempt  was  made  to  conform  to  ASA  Standard  FORTRAN  [6]  as  far  as 
possible.  Specific  areas  where  the  code  is  known  to  depart  from  Standard  FORTRAN  are: 

(1)  The  use  of  subscripted  integer  variables  as  subscripts 

(2)  The  ‘OPEN'  instruction,  which  sets  up  disc  files  for  input  and  output 

(3)  The  'TYPE'  statement,  which  is  used  to  output  diagnostic  information  directly  to  a 
terminal. 

In  addition,  certain  tests  to  allow  for  rounding  errors  involve  constants  appropriate  to  a 36-bit 
processor. 

The  plotting  program,  SPP,  on  the  other  hand,  is  rather  more  dependent  on  the  software 
available  at  this  particular  installation.  In  particular,  it  makes  use  of  system  subroutines  which 
drive  a digital  incremental  plotter. 

On  the  DECsystem-10  with  a Kl  processor,  run  times  are  about  2 sec.  to  calculate  the 
coordinates  for  a 200-point  projection,  and  of  the  order  of  10  sec.  to  create  the  plot  file. 

Program  listings  have  been  deposited  with  the  ARL  Computer  Centre  and  are  available 
within  the  ARL  system  from  DECtape  D/933. 

3.  THEORY 

For  a lattice  with  parameters  a,b,c.*.f$,y.  the  angle  <b  between  directions  wrir  and  i/rV  is  [7] 

cos  <f>  = {uu'a2  + vv'b2  + ww'c2  + (v\\  ' + wv')bc  cos  * + (ie u'  + uw')ca  cos  /S 

+ (uv'  + vu')ab  cos  y]/DD'  (I) 
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where 

D = (u2a2  + v2b2  + tv2c2  + 2 vwbc  cos  a + 2iv«ca  cos  fi  + luvab  cos  y)*  (2) 

is  the  repeat  distance  along  the  direction  uvw , and  D',  correspondingly  defined,  the  repeat  distance 
along  the  direction  u'v'w'. 

Because  of  the  relationship  between  directions  in  the  rer! lattice  and  planes  in  the  reciprocal 
lattice,  and  vice-versa,  4>  is  also  the  angle  between  the  normals  to  lattice  planes  uvw  and  u'v'w'  in 
the  reciprocal  lattice.  Conversely,  by  replacing  the  parameters  a.b.c.a.fi.y  by  the  reciprocal  lattice 
parameters  a*,b*,c*,ai*,f)*,y*,  <f>  becomes  the  angle  between  the  normals  to  the  planes  uvw  and 
u'v'w'  in  the  real  lattice. 

The  reciprocal  lattice  parameters  are  defined  by 

a*  = (he  sin  x)jV  cos  a*  = (cos  fi  cos  y — cos  a)/(sin  /3  sin  y) 

b*  = ( ca  sin  fi)/V  cos  /3*  = (cos  y cos  a-cos  y3)/(sin  y sin  *) 

c*  = (ab  sin  y)/V  cos  y*  = (cos  a cos  /?  — cos  y)/( sin  a sin  /?) 

where 

V = abc  ( I — cos  2a  — cos  2|8  - cos  2y  + 2 cos  a cos  cos  y)  * 
is  the  volume  of  the  unit  cell  [8]. 

From  equation  (2),  it  can  be  shown  to  be  necessary  (but  not  sufficient)  that 

(1)  For  D to  be  a maximum, 

h-  = 0 (3a) 

(2)  For  »■  fixed,  D will  be  a maximum  for 

v = (m  c//>)(cos  (3  cos  y — cos  a)/(sin2y)  (3b) 

(3)  For  **■  and  v fixed,  D will  be  a maximum  for 

ii  = (rb  cos  y-t-ii  c cos  f))ja  (3c) 

Consider  a reference  sphere  centred  at  the  origin  of  the  lattice  (Figure  la).  A point  Q on  the 
surface  of  the  sphere  uniquely  represents  a lattice  direction  passing  through  the  origin  and  is 
represented  stereographically  on  the  horizontal  plane  through  the  centre  of  the  sphere  by  the 
point  P.  which  is  the  point  where  this  plane  is  intersected  by  the  line  connecting  Q to  the  south 
pole.  This  plane  is  the  projection  plane  (Figure  lb). 

Consider  three  reference  directions  (31,  Q 2 and  Q 3.  defined  by  Miller  indices  wl.rl.nl, 
u2.r2.w2  and  u3.r3.w3,  respectively.  Q 1 defines  the  north  pole  and  becomes  the  central  pole  on 
the  projection.  Q2  and  £33  are  perpendicular  to  Q I and  to  each  other  and  therefore  lie  on  the 
perimeter  of  the  projection.  We  may  define  a coordinate  system  R.0.<f>  using  these  directions  as 
coordinate  axes:  Q 1 (<f>  = 0 ),  Q2  (<t>  = 90’,  0 — 0°),  (33  (<f>  = 90’,  6 = 90n).  The  arbitrary 
direction  Q.  u.v.w.  makes  angles  of  <f>.  if  and  y with  £31,  Q2  and  (33.  These  angles  may  be  deter- 
mined using  equation  (1).  The  coordinate  0 in  the  system  defined  above  is  then  given  by 

cos  0 = (cos  i/i)/( sin  <f> ) 

the  positive  or  negative  value  of  0 being  taken  depending  on  whether  y is  greater  or  less  than  90  . 

The  coordinate  r is  the  distance  of  the  point  P from  the  centre  of  the  projection  and  has  the 
value  R tan  612.  where  R is  the  radius  of  the  reference  sphere. 

4.  PROGRAM  SPC— STEREOGRAPHIC  PROJECTION  COORDINATES 

The  program  is  written  as  a collection  of  subroutines,  which  are  called  by  the  main  program 
in  the  appropriate  order.  Extensive  use  is  made  of  COMMON  to  minimise  data  transfers  and 
core  storage.  The  program  allows  data  to  be  determined  for  up  to  222  points.  This  number  of 
points  has  been  found  ample  for  most  applications.  However,  the  combination  of  digits  '222'  is 
unique  to  this  context  and  the  program  can  be  edited  to  handle  any  desired  maximum  by  globally 
changing  '222'  to  the  desired  figure  before  compiling. 

The  program  basically  considers  crystal  directions,  but  an  appropriate  request  in  the  input 
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data  causes  the  lattice  to  be  replaced  by  the  corresponding  reciprocal  lattice.  Operation  of  the 
program  on  the  reciprocal  lattice  then  gives  data  which  relate  to  lattice  plane  normals  in  the 
direct  lattice.  In  this  documentation  of  the  program,  all  reference  is  to  directions.  A corresponding 
description  applies,  however,  for  lattice  plane  normals.  It  is  not  possible  to  consider  both  planes 
and  directions  simultaneously. 

4.1.  Description  of  the  Program 

4.1.1.  Main  Program 

Subroutines  are  called  in  sequence. 

4.1.2.  Subroutine  INPUT 

Disc  files  are  opened  for  input  and  output.  Input  data  arc  read  from  file  and  their  plausibility 
checked.  The  main  data  required  are: 

(1)  Title 

(2)  Lattice  parameters 

(3)  Information  defining  the  centre  of  the  projection 

(4)  Information  defining  the  range  of  points  to  be  accepted 

(5)  Whether  planes  or  directions  are  required 

(6)  Radius  of  projection 

(7)  Maximum  number  of  points  allowed. 

The  input  data  are  also  used  to  specify  whether  stereographic,  polar,  or  cartesian  coordinates  are 
required  for  the  output.  Detailed  specifications  for  the  data  required  are  included  as  comments 
in  the  subroutine. 

4.1.3.  Subroutine  RECLAT 

The  reciprocal  lattice  is  set  up  corresponding  to  the  direct  lattice  specified.  The  direct  lattice 
is  replaced  by  the  reciprocal  lattice  if  lattice  plane  normals  are  to  be  considered.  Constants  to  be 
used  in  other  subroutines  are  determined. 

4.1.4.  Subroutine  ORTHOG 

From  the  input  data  defining  the  orientation  of  the  projection,  the  program  sets  up  reference 
directions  at  the  centre  of  the  projection  (01).  and  on  the  circumference  at  0 = 0 and  90° 
(02.03,  respectively).  The  following  cases  are  distinguished: 

If  only  one  direction  is  given,  it  is  taken  as  01.  If  two  directions  are  given  and  they  are 
perpendicular  to  one  another,  they  are  taken  as  01  and  02.  respectively.  If  two  directions  are 
given  and  they  are  not  perpendicular,  01  is  set  up  perpendicular  to  them  both.  i.e.  these  two 
directions  are  regarded  as  lying  on  the  perimeter  of  the  projection. 

4.1.5.  Subroutine  RATNL 

An  attempt  is  made  to  rationalise  indices  of  directions  set  up  within  subroutine  ORTHOG. 

Additional  constants  to  be  used  in  other  subroutines  are  determined. 


4.1.6.  Subroutine  LIMITS 

It  was  originally  envisaged  that  ‘low-index’  directions  could  usefully  be  displayed  more 
prominently  on  a projection,  since  these  are  usually  of  greater  importance.  In  generalising  to 
non-cubic  lattices,  this  criterion  was  no  longer  satisfactory  and  a new  rule-of-thumb  was  adopted, 
viz.  that  the  smaller  the  repeat  distance.  D,  along  a given  direction,  the  more  ‘important'  it  was 
likely  to  be.  Thus,  it  was  decided  to  classify  the  points  obtained  into  groups  depending  on  their 
D-para  meters. 

A series  of  n sets  of  Miller  indices  may  be  specified  in  the  input  (I  ^ n sg  10).  The  repeat 
distance  Dn  is  calculated  for  each  set.  A candidate  pole  is  then  included  in  the  output  if  its  repeat 


distance  D lies  within  the  range  D 1 § D ^ Dn.  Acceptable  values  of  D,  moreover,  may  be 
separated  into  n-  1 groups,  depending  on  the  value  of  D. 

This  structure  facilitates  obtaining  a projection  with  a convenient  number  of  points,  especially 
for  low-symmetry  lattices.  If  the  range  of  calculation  is  defined  as  a multiple  range  (i.e.  n > 2), 
and  the  number  of  points  requested  is  exceeded  during  computation,  the  program  deletes  the 
uppermost  range  and  starts  again.  Since  non-integral  indices  may  be  used  here,  it  is  possible  to 
define  a very  narrow  range,  so  that  only  a single  type  of  direction  is  included.  This  facility  is  used 
for  preparing  ideal-texture  pole  figures  (Section  5.12). 

4.1.7.  Subroutine  GENUVW 

This  subroutine  constitutes  the  core  of  the  program.  Sets  of  indices  t/,r,tr  are  generated, 
tested  to  see  whether  they  fall  within  the  desired  constraints,  and  the  coordinates  of  the  direction 
represented  are  calculated  in  relation  to  the  three  reference  directions.  An  algorithm  was  required 
for  generating  suitable  sets  of  indices  w.r.tr.  For  orthogonal  lattices,  simple  stepping  of  u , v and  tr 
on  either  side  of  zero  would  be  efficient,  but  for  crystals  of  lower  symmetry,  various  problems 
arise  and  the  following  procedure  was  developed. 

As  discussed  above,  the  repeat  distance  is  used  as  the  criterion  for  inclusion  in  the  output. 
For  given  values  of  r and  >r.  there  will  be  some  value  I/O  which  will  lead  to  a minimum  of  D. 
This  can  be  determined  by  differentiating  the  expression  for  D (Equation  3c).  For  these  values  of 
v and  tr,  values  of  u can  then  be  stepped  out  on  either  side  of  uO.  the  procedure  terminating  when 
the  combinations  of  u.  r and  it  cease  to  be  acceptable.  A similar  process  occurs  for  each  com- 
bination of  c and  tr.  A related  procedure  is  used  for  each  tr  to  determine  a value  iO  from  which 
values  of  r are  stepped  (Equation  3b).  Values  of  it  are  simply  stepped  about  zero  (Equation  3a). 

The  coordinates  of  the  point  i/.r. tr  are  determined  as  0 and,  depending  on  the  input 
parameters,  r or  <f>. 

As  acceptable  points  are  found,  the  related  data  are  stored.  Because  values  of  0 greater  than 
180  are  specially  considered  in  the  program,  the  results  are  conveniently  stored  in  two  blocks, 
thereby  reducing  the  sorting  required  later.  If  the  available  storage  becomes  full,  execution 
'erminates  unless  a multiple  range  has  been  defined  for  the  range  of  calculation  (see  Section 
4.1.6),  in  which  case,  the  least  significant  range  is  discarded  and  the  calculation  repeated  with  the 
reduced  range. 

4.1.8.  Subroutine  ORDER 

The  stored  data  are  sorted  into  increasing  order  with  respect  to  0,  and.  for  each  f>.  in  de- 
creasing order  with  respect  to  r or  <b. 

4.1.9.  Subroutine  HEAD 

General  information  is  presented  on  a title  page. 

4.1.10.  Subroutine  OUTPUT 

If  required,  all  relevant  data  are  dumped  on  a binary  file.  Hard-copy  output  is  initially  sent 
to  file  in  a form  determined  by  the  input  parameters. 

4.2.  Operation  of  the  Program 

Detailed  specifications  for  the  data  required  for  program  SPC  are  given  as  comments  in 
Subroutine  INPUT.  Table  1 gives  a sample  set  of  input  data.  It  requests  a 111  projection  of 
lattice  planes  for  a hexagonal  lattice,  eja  = 1 - 59.  All  planes  are  to  be  included  whose  interplanar 
spacing  is  greater  than  that  of  the  400  planes.  The  output  will  give  stereographic  coordinates  <t>, 
ft,  and  a binary  file  will  be  produced  for  the  plotting  program. 

The  teletype  record  when  running  the  program  is  shown  in  Table  2.  The  number  of  points 
specified  (222  by  default)  was  exceeded.  However,  since  a multiple  range  was  specified,  the  least 
significant  range  (300  400)  was  discarded  and  computation  repeated. 

Part  of  the  output  is  shown  in  Table  3. 
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5.  PROGRAM  SPP — Stercographic  Projection  Plotter 


5.1.  Description  of  the  Program 


5.1.1.  Standard  Projections 


This  program  accepts  as  input  the  coordinate  data  produced  by  program  SPC  and  plots  the 
projection  on  a digital  incremental  plotter.  Other  parameters  are  input  separately  to  determine 
the  detailed  format. 

The  program  makes  extensive  use  of  a suite  of  system  subroutines  which  draw  lines  and 
symbols  on  the  plotter.  Such  subroutines  are  available  at  most  installations  which  use  graph- 
plotting facilities.  However,  the  coding  in  SPP  may  need  to  be  modified  before  the  progcom  can 
be  used  on  some  other  system.  ./ 

A specification  of  the  operations  performed  by  these  subroutines  is  given  in  Appendix. 
The  prc.ram  permits  standard  projections  up  to  21  inches  in  diametpri  although  for 
diameters  greater  than  1 1 inches,  or  if  requested  specifically,  the  projectkxf  is  plotted  as  two 
semi-circles.  These  restrictions  are  necessitated  by  the  limited  plotter  whfitff  If  requested,  only  the 
first  quadrant  will  be  plotted.  / ' 

For  hexagonal  lattices,  it  is  more  convenient  to  use  the  Mitfer-Bravais  4-index  notation.  The 
conversion  is  executed  as  the  default  option,  but  may  Nxfippressed  if  desired. 


5.1.2.  Ideal  Sheet  Textures 


An  additional  facility  is  the' plotting  of  pole  figures  of  so-called  "ideal"  sheet  textures  of 
rolled,  or  rolled  and  recrystallised,  metal  sheets  [9],  Because  of  the  symmetry  of  the  rolling 
process,  the  pole  figure  of  such  a sheet  is  symmetrical  about  the  rolling  direction.  For  reversed- 
rolling  at  least,  it  is  also  symmetrical  about  the  transverse  direction.  Accordingly,  an  ideal 
texture  can  be  simulated  by  reflecting  the  pole  figure  through  the  radii  representing  the  rolling 
and  transverse  directions. 

The  program  arranges  for  all  points  in  the  pole  figure  to  be  reflected  into  the  first  quadrant. 
By  adding  a rotation  to  the  source  data  beforehand,  the  rolling  direction  can  be  varied.  The 
program  produces  a set  of  quadrant  pole  figures  for  a specified  set  of  rotations. 

The  procedure,  although  somewhat  naive,  is  simple  in  concept  and  can  be  of  great  use  in 
identifying  texture  components. 


5.2.  Operation  of  the  Program 

5.2.1.  Data  Required 

File  SPC. BIN  binary  output  file  from  program  SPC. 

File  SPP. DAT  auxiliary  data  defining  the  format  of  the  projection: 

( 1 ) Diameter  of  reference  circle 

(2)  A code  defining  the  size  indices  arc  to  be  plotted 

(3)  A code  requesting  that  a circle,  two  semi-circles,  or  a quadrant  be  drawn 

(4)  A code  limiting  the  number  of  data  points  to  be  plotted 

(5)  A code  defining  the  angle  at  which  the  indices  are  to  be  plotted 

(6)  A code  suppressing  the  conversion  to  Millcr-Bravais  indices 

(7)  Three  parameters  defining  a sequence  of  angles  for  generating  ideal-texture  pole  figures. 

Detailed  descriptions  of  the  data  required,  and  their  format,  are  included  as  comments  in  the 
program. 

5.2.2.  Examples 

The  results  from  the  earlier  test  (shown  in  Table  3)  were  plotted  with  all  zeros  for  the 
auxiliary  data  (thereby  implementing  default  settings  throughout).  The  resulting  plot  is  shown 
in  Figure  2. 


As  an  example  of  the  meal-texture  facility,  the  5th  line  of  the  data  (Table  1)  was  changed 
to  read 

0-99  00  I 01  0 0 

i.e.  so  as  to  restrict  the  output  to  include  only  1 0 0 planes.  After  processing  with  programs 
SPC  and  SPP.  and  with  the  auxiliary  data 

000000  150  80 

the  set  of  ideal  textures  show  n in  Figure  3 was  obtained.  (Note  that  when  the  ideal-texture  option 
is  being  used,  many  of  the  auxiliary  data  items  need  not  be  set  because  they  are  reset  to  standard 
values  by  the  program.) 

6.  CONCLUSION 

Two  FORTRAN  programs  have  been  written  for  generating  the  points  on  standard  stereo- 
graphic  projections  for  the  most  general  of  cases  and  for  plotting  these  projections.  The  programs 
have  operated  reliably  over  several  years  and  have  been  used  for  drawing  projections  for  use  in 
studies  on  deformation,  texture,  and  electron  microscopy  and  diffraction.  Idealised  pole  figures 
of  sheet  textures  may  also  be  prepared. 
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APPENDIX 


Specification  of  subroutines  used  to  drive  the  CALCOMP  digital  incremental  plotter.  Only 
those  facilities  are  described  which  are  invoked  by  the  program  SPP.  The  specification  has  been 
adapted  from  the  complete  specification  of  the  set  of  plotting  subroutines  ALLPLT  [10], 

Subroutine  PLOT 

Purpose:  To  provide  the  basic  movements  of  the  plotter  pen. 

Usage:  CALL  PLOT(LU,X,Y,N) 

Description  of  Parameters: 

LU:  Integer  input.  Logical  unit  number. 

X,Y:  Real  inputs.  X and  Y coordinates  of  the  point  to  be  plotted. 

N:  Integer  input. 

N = 1:  X.Y  are  scaling  factors  (default:  1 .,!.)• 

N = 2:  Reset  origin.  The  current  position  of  the  pen  is  defined  as  (X,Y). 

N = 3:  Pen  is  lifted,  moved  to  X.Y  and  dropped  (not  used). 

N = 4:  Pen  is  dropped  and  moved  to  X.Y. 

N = 5:  Pen  is  lifted,  moved  to  X.Y,  but  not  dropped. 

Remarks:  The  positive  X-direction  is  along  the  paper  towards  the  feed  roll.  The  positive  Y- 
direction  is  from  right  to  left  parallel  to  the  drum  axis  when  looking  at  the  plotter.  The  initial 
position  of  the  pen  is  at  the  centre  of  the  paper  and  this  is  initially  defined  as  the  origin. 

Subroutine  WHERE 

Purpose:  To  return  the  current  position  of  the  pen  relative  to  the  current  origin  (in  units  con- 
sistent with  the  current  scale  factors). 

Usage:  CALL  WHERE  (LU.X.Y). 

Description  of  Parameters: 

LU:  Integer  input.  Logical  unit  number. 

X.Y : Real  output.  X and  Y coordinates  of  current  pen  position. 

Subroutine  SYMBOL 

Purpose:  To  draw  characters  on  the  plotter. 

Usage:  CALL  S YM BOL(LU. X.Y. HEIGHT, ASCI I.TH ETA. N) 

Description  of  Parameters: 

LU:  Integer  input.  Logical  unit  number. 

X.Y.  (1)  N > 0:  coordinates  of  the  lower  left  hand  corner  of  the  character  string  to  be 
plotted.  On  exit,  the  pen  is  left  at  the  lower  left-hand  corner  of  the  next  character  to  be 
plotted  (pen  up). 

(2)  N < 0:  Coordinates  of  centre  of  special  symbol.  On  exit,  the  pen  is  left  at  X.Y  (pen 
up). 

X and  Y may  be  real  or  integer,  indicating,  respectively,  coordinates  in  inches  with  respect 
to  the  current  origin,  or  coordinates  in  0 01  inch  increments  with  respect  to  the  current 
position  of  the  pen. 

HEIGHT:  Real  input.  The  height  in  inches  of  the  characters  to  be  drawn. 

ASCII.  (I)  N > 0:  ASCII  is  the  first  word  of  the  string  of  characters  to  be  plotted. 

(2)  N < 0:  ASCII  must  be  an  integer.  The  only  values  called  here  are  N = — 1, 
ASCII  = 0,  which  request  a square  with  a diagonal  cross. 

THETA:  Real  input.  Angle  (in  degrees)  the  string  of  characters  is  to  make  with  the  X-axis. 
THETA  increases  in  the  anticlockwise  direction. 
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N:  (1)  N > 0:  N characters  are  plotted. 

(2)  N < 0:  The  special  symbol  specified  in  ‘ASCII’  is  to  be  drawn. 

Subroutine  NUMBER 

Purpose:  To  draw  the  value  of  a real  number  on  the  plotter. 

Usage.  Call  NUMBER(LU,X,Y,HEIGHT,FLOAT,THETA,N) 

Description  of  Parameters: 

LU,X,Y,HEIGHT,THETA:  As  for  subroutine  SYMBOL. 

FLOAT : the  real  number  whose  value  is  to  be  drawn. 

N:  Number  of  decimal  places  to  be  retained.  The  special  case  N = -1  prints  only  the 
integral  part  of  N,  without  a decimal  point. 

Remarks:  Leading  zeros  are  omitted  and  the  number  is  left-justified.  If  FLOAT  is  less  than  0-5 
a single  zero  is  drawn.  ' ’ 
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TABLE  I 

Sample  data  for  program  SPC 

ALPHA  TITANIUM 
2950  2-950  4-683 

90-  90-  120 

1 1 1-1  10 
000200300 
PLANES 


I 


TABLE  2 

Teletype  record  for  running  program  SPC 

RUN  SPC 

STORAGE  EXHAUSTED 

restarted  with  reduced  range 

104  POINTS 
end  of  execution 

CPU  time:  I 40  elapsed  time:  2-78 
exit 


j 


TABLE  3 

Hard-copy  output  from  program  SPC  using  data  in  Table  1 

(The  layout  has  been  somewhat  condensed) 


ALPHA  TITANIUM 

DIRECT  LATTICE. 

A =2-950 
ALPHA  = 90  000 

RECIPROCAL  LATTICE: 
A*  =0-391 
ALPHA*  90  000 


B 2-950 
BETA  =90-000 


B*  =0-391 
BETA*  90-000 


C 4-683 

GAMMA  = 120-000 


C*  = 0-214 
GAMMA*  = 60-000 


COORDINATE  DATA  FOR  STANDARD  (1,1,1)  PROJECTION  FOR  LATTICE  PLANE 
NORMALS. 

RANGE  OF  D-PARAMETER:  0.0  0.0  0.0  3.0  0.0  0.0  104  POINTS. 

STEREOGRAPHIC  COORDINATES  THETA,  PHI. 

CENTRE  POLE:  1-00  1 -00  1-00  THETA  0:  -I  00  1-00  0 00 


THETA 

U 

V 

W 

PHI 

U 

V 

w 

PHI 

U 

V 

0-00 

-1 

1 

0 

90 

00 

-1 

3 

1 

47-76 

0 

2 

5-94 

-2 

3 

I 

68 

50 

9-84 

-1 

2 

1 

56 

97 

14-58 

-2 

2 

1 

85 

47 

-1 

3 

2 

46-23 

17-34 

-2 

3 

2 

66 

22 

27-49 

-1 

1 

1 

81 

73 

-1 

2 

2 

55-69 

0 

1 

37-97 

-2 

2 

3 

79 

03 

40-93 

-l 

2 

3 

56 

43 

46-14 

-1 

1 

2 

77 

20 

0 

2 

3 

33-95 

50-52 

-I 

2 

4 

57 

87 

57-36 

-1 

1 

3 

75 

15 

0 

1 

2 

38  ■ 78 

64-34 

-1 

1 

4 

74 

15 

68-98 

-1 

1 

5 

73 

62 

0 

1 

3 

46-60 

74-65 

0 

1 

4 

51 

89 

77-95 

0 

1 

5 

55 

50 

90  00 

0 

0 

1 

72 

52 

1 

1 

4 

34-08 

1 

1 

102-05 

1 

0 

5 

55 

50 

105-35 

1 

0 

4 

51 

89 

111-02 

1 

-1 

5 

73 

62 

1 

0 

3 

46-60 

115  66 

l 

-1 

4 

74 

15 

122-64 

1 

-1 

3 

75 

15 

1 

0 

2 

38-78 

129-48 

2 

-1 

4 

57 

87 

133-86 

1 

-1 

2 

77 

20 

2 

0 

3 

33-95 

139-07 

2 

-1 

3 

56 

43 

142-03 

2 

-2 

3 

79 

03 

152-51 

1 

-1 

1 

81 

73 

2 

-1 

2 

55-69 

1 

0 

162-66 

3 

-2 

2 

66 

22 

165-42 

2 

-2 

1 

85 

47 

3 

-1 

2 

46-23 

170-16 

2 

-1 

1 

56 

97 

174-06 

3 

-2 

1 

68 

50 

W 

1 


PHI 

28-84 


U 

1 


V W 
1 1 


PHI 

0 00 


U V W PHI 


1 29  66  1 2 2 11-36 


3 25-89  1 


2 14-73 


1 29  66  2 1 2 11-36 


Table  3 — continued 


180-00 

1 

-1 

0 

90-00 

3 

-1 

1 

185-94 

3 

-2 

0 

71-81 

189-84 

2 

-1 

0 

61-52 

194-58 

3 

-1 

0 

51  36 

197-34 

3 

-2 

-1 

75-81 

207-49 

3 

-2 

-2 

79  97 

2 

-1 

-1 

217-97 

3 

-1 

-2 

62-44 

220-93 

2 

-1 

-2 

75  91 

226  14 

2 

0 

-1 

44  12 

230-52 

2 

-1 

-3 

81-93 

237-36 

2 

-l 

-4 

86-49 

1 

0 

-1 

244-34 

2 

0 

-3 

63-32 

248-98 

1 

0 

-2 

70  33 

2 

1 

-1 

254-65 

1 

0 

-3 

79-96 

2 

1 

-2 

257-95 

1 

0 

-4 

85  93 

260-09 

1 

0 

-5 

89-87 

270  00 

1 

1 

-4 

69  04 

1 

1 

-3 

279  91 

0 

1 

-5 

89-87 

282-05 

0 

1 

-4 

85  93 

285-35 

0 

1 

-3 

79-96 

I 

2 

-2 

291 -02 

0 

1 

-2 

70-33 

1 

2 

-1 

295-66 

0 

2 

-3 

63-32 

302-64 

-1 

2 

-4 

86  49 

0 

1 

-1 

309-48 

-1 

2 

-3 

81  93 

313-86 

0 

2 

-1 

44- 12 

319  07 

-1 

2 

-2 

75-91 

322-03 

-1 

3 

-2 

62-44 

332-51 

_2 

3 

-2 

79-97 

-1 

2 

-1 

342  66 

-2 

3 

-1 

75-81 

345-42 

-1 

3 

0 

51-36 

350  16 

-1 

2 

0 

61-52 

354-06 

-2 

3 

0 

71-81 

47-76  2 0 I 28-84 

68-62  3 -I  -I  56-57  1 0 0 34-31 

54-46  2 l 0 20-51 

31  06 
41  29 

60-86  1 1 -2  49-69  1 1 -1  34-97 

41-29 
31  -06 

54-46  I 2 0 20-51 

68-62  -I  3 -I  56-57  0 1 0 34-31 


2 1 


1 I 


1 2 


1 12-04 


0 17-48 


1 12-04 


FIGURE  2.  STANDARD  1121  PROJECTION  FOR  LATTICE  PLANES  IN  TITANIUM. 


FIGURE  3.  SIMULATED  100  POLE  FIGURE  FOR  TITANIUM  LATTICE.  SEVERAL 
TEXTURES  OF  THE  TYPE  (111)  [UVW]  ARE  REPRESENTED. 
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